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FOREWORD 


The  research  described  in  this  report,  A  Survey  of  Techniques  for 
the  Analysis  of  Sampled-Data  Systems  with  a  Variable  Sampling  Rate,  by 
George  A.  Bekey,  was  carried  out  under  the  technical  direction  of 
C.  T.  Leondes  and  G.  Estrin  and  is  part  of  the  continuing  program  in 
Control  Systems  Theory. 

This  project  is  conducted  under  the  sponsorship  of  Aeronautical 
Systems  Division,  Flight  Control  Laboratory,  Project  Engineer 
Charles  Harmon,  Wright-Patterson  Air  Force  Base.  Submitted  in  partial 
fulfillment  of  Contract  Number  AF  33(616 ) -71 39,  Project  Number  8225. 
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ABSTRACT 


» 

This  report  presents  several  techniques  which  can  be  used  for  the 
analysis  of  sampled-data  systems  with  a  non-constant  sampling  period.  It 
is  shown  that  the  application  of  z-transform  techniques  is  limited  to  cases 
where  the  sampling  pattern  repeats  periodically.  Several  special  cases, 
including  the  cyclically-varying  sampling  period  and  the  "Skipped  sample" 
problem  are  outlined.  The  most  general  methods  available  are  based  on  a 
direct  solution  of  the  system  difference  equations.  These  equations,  while 
time  consuming  to  solve,  do  make  possible  the  evaluation  of  transient  re¬ 
sponse  sample-by-sample.  The  advantages,  limitations  and  possible  exten¬ 
sions  of  the  various  methods  are  outlined.  The  report  includes  a  number  of 
simple  examples  and  an  extensive  bibliography. 
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1. 


INTRODUCTION 


The  analysis  of  sampled-data  systems  with  a  fixed  sampling  period  has 
been  treated  in  a  large  number  of  publications,  primarily  by  use  of  the 

z-transform  method.  Only  a  relatively  small  amount  of  attention  has  been 
devoted  to  the  larger  class  of  sampling  systems  in  which  the  sampling  inter¬ 
val  is  variable,  depending  either  on  time  or  on  a  system  variable;  the  former 
being  linear  and  the  latter  non-linear.  Sampled-data  systems  with  a  variable 
sampling  rate  show  promise  of  having  two  extremely  interesting  properties: 

(1)  The  operation  of  a  system  at  its  minimum  sampling  rate,  consis¬ 
tent  with  required  performance  and  stability  criteria,  leads  to  a 
class  of  adaptive  systems  in  which  the  energy  or  control  effort  is 
minimized.  Such  control  systems  may  be  extremely  important 
for  applications  when  available  energy  is  severely  limited,  such 
as  satellite  attitude  control. 

(2)  The  operation  of  a  system  which  adjusts  its  sampling  rate  accord¬ 
ing  to  performance  requirements  may  closely  resemble  the  behav¬ 
ior  of  a  human  operator  and  therefore  may  lead  to  clues  about 
certain  functions  of  the  central  nervous  system. 

The  purpose  of  this  paper  is  to  review  and  summarize  the  major  tech¬ 
niques  described  in  the  literature  for  the  analysis  of  sampled-data  systems 
with  a  non-constant  sampling  interval.  A  general  and  unified  approach  to  this 
problem  seems  to  be  lacking.  Most  of  the  available  papers  concerning  vari¬ 
able  sampling  assume  either  periodic  variations  in  sampling  rate  or  random 
variations  in  sampling  intervals.  This  report  presents  an  outline  of  the  major 
methods  presently  available  and  discusses  briefly  their  applicability  and  limi¬ 
tations. 

A  recent  article  by  Jury  presents  a  brief  discussion  of  various  types  of 
22 

sampling  schemes. 

1,  1  Classification 

Sampled  data  systems  with  a  non-constant  sampling  rate  may  be  divided 
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into  three  general  types: 


(1)  Time-dependent  sampling:  The  sampling  intervals  are  given  by 
a  known  function  of  time.  Nearly  all  the  literature  on  variable¬ 
sampling  rate  systems  is  concerned  with  this  type  of  sampling. 

Specific  systems  include:  (a)  Cyclic-rate  systems  (where  the 

sampling  interval  is  a  periodic  function  of  time,  i.  e. ,  a  pattern  ! 

of  sampling  intervals  repeats  periodically);  (b)  Skip-sampled 

systems  (in  which  samples  in  an  otherwise  constant  rate  system 

are  occasionally  missed);  and  (c)  Piecewise-constant  sampling 

rate  systems  (where  the  sampling  rate  changes  discontinuously 

after  a  number  of  samples,  remaining  constant  between  changes). 

(2)  Nonlinear  sampling:  The  sampling  intervals  are  functions  of  de¬ 
pendent  variables  or  system  performance.  Clearly,  this  prob¬ 
lem  is  considerably  more  difficult  than  that  of  time  dependent  sam¬ 
pling,  and  very  little  information  is  available  in  the  literature  on 
methods  of  handling  such  problems. 

(3)  Random  sampling:  In  this  type  of  system  the  sampling  interval  is 
a  random  variable.  The  synthesis  of  optimum  filters  for  such 

systems  with  stationary  random  sampling  has  been  treated  in  the 

,  1,11,15 

literature. 

In  the  present  report  we  shall  concern  ourselves  exclusively  with  deter¬ 
ministic  sampling  and  omit  any  further  discussion  of  random  sampling.  Further¬ 
more,  the  report  is  concerned  with  constant  pulse  width  systems,  preferably 
those  which  can  be  analyzed  using  the  impulse  approximation.  A  number  of 
extremely  interesting  papers  have  dealt  with  pulse-width  modulation  in  sam- 
pled-data  systems  but  this  problem  will  be  considered  as  being  outside  the  » 

scope  of  this  discussion. 

1.  2  Methods  of  Analysis 

A  number  of  methods  have  been  proposed  for  the  analysis  of  sampled- 
data  systems  with  variable  sampling  rate.  In  many  cases,  the  methods  are 
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applicable  only  to  the  case  of  periodic,  time  dependent  sampling  rates.  The 
major  techniques  known  to  the  author  fall  into  the  following  classification: 

(1)  The  solution  of  difference  equations  with  variable  coefficients  in 
the  time  domain. 

(2)  The  use  of  z-transform  and  modified  z-transforms. 

(3)  The  use  of  matrix  methods  for  formulation  of  the  solution. 

(4)  The  application  of  state- variable  concepts  for  a  general  formula¬ 
tion  of  the  problem. 

(5)  Approximate  solutions,  using  continuous  system  techniques  to 
approximate  sampled  data  systems. 

The  five  methods  above  will  be  examined  in  small  detail  in  the  following 
paragraphs.  While  it  is  clear  from  the  discussion  which  follows  that  there  is 
a  great  deal  of  overlap  between  these  five  methods,  they  are  presented  sepa¬ 
rately  for  convenience  and  because  they  appear  separately  in  the  literature. 

2.  ANALYSIS  USING  DIFFERENCE  EQUATIONS  WITH  TIME -VARYING 

COEFFICIENTS 

The  time-domain  behavior  of  a  sampled-data  system  can  be  described 
using  difference  equations,  even  if  the  sampling  intervals  are  not  of  constant 
length.  Thus  a  q-th  order  system  can  be  described  by  q  1-st  order  difference 
equations  or,  by  substitution,  a  single  q-th  order  equation  which  expresses 
the  output  at  any  sampling  instant  in  terms  of  the  input  and  output  values  at 
the  (q-1)  past  sampling  instants. 

Consider  the  error-sampled  system  of  Figure  1,  which  includes  a  zero- 
order  hold  and  a  linear  "plant”  described  by  its  transfer  function  G(s). 


ERROR-SAMPLED  SYSTEM 
FIGURE  1 
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The  sampler  is  assumed  to  close  af  instants  t^^,  ,  where  the  sampling 

intervals 

T  =  t  ^  -  t 
n  n+1  n 

are  not  necessarily  equal.  For  the  moment,  the  sampling  instants  will  be 
assumed  known,  and  the  removal  of  this  restriction  will  be  discussed  later. 
The  output  of  the  system  during  the  interval 


t  <  t  <  t 
n  —  —  n+1 


is  then  given  by  the  relation® 


c(t)  =  X^g(t-t^)  + 


q-1 

y  f  (t-t  ) 

Zy  P  n 


p=0 


d'P'cCt) 

1 

dt<p’ 

II 

’h* 

n 


(la) 


where  t^  =  time  at  the  start  of  the  (n+1)  sampling  interval  of  duration  T^ 

X  =  output  of  zero-order  hold  =  e(t  ) 
c(t)  =  continuous  system  output 
g(t)  =  step  function  response  of  system  G(s) 
q  =  order  of  the  system 

and  the  f  (t)  are  time  functions  which  result  from  the  initial  conditions  at  the 
P 

start  of  the  n-th  sampling  interval. 


Since  the  system  is  of  order  q,  we  need  q  first  order  equations  to  com¬ 
pletely  describe  it.  The  additional  equations  can  be  obtained  by  successive 
differentiations  of  Equation  (la)  to  give: 


dc 

dt 


t  <t<t  , 
A - n+1 


=  X  g^^^t-t  )  + 
n  n 


(q-1) 

I 


p=0 


- 1 

a 

1 

a 

t=t  - 

n 


(lb) 


dl^ 


t  <t<t  , 
n - n+1 


(q-i) 

I' 

p-0 


P  n 


d<P>c(t) 

dt<p> 

t=t  . 

(Ic) 
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if  we  define 


dt 


(i) 


A  c(i) 

n 


t=t 


n 


and  let  the  q  equations  reduce  to  the  following  set  of  difference 

equations: 


q-1 


-C  )g  (T  )  +  y  UT 
n+1  nn  n  ^pnn 


(P) 


p=0 


(2a) 


q-1 

C"‘r'=  (R„-C^)g'^ 


p  n  n 


(2q) 


p=0 


where  we  have  let  -  C^.  Then  the  system  is  completely  described  by 

Equations  (2a)  to  (2q). 

The  above  equations  can  be  replaced  by  an  equivalent  set  which  does  not 
require  knowledge  of  the  derivatives  of  the  output  at  the  sampling  instants  but 
only  the  value  of  the  output.  The  equations  are  written  in  the  form: 

q-1 


'n+1 


=  (R  “C  )g(T  )  +  y  f  (T  )C^P^ 
nn'’n  Z^pnn 


p=0 


(3a) 


q-1 


n 


(R 


,-C  Jg(T  .)  +  y  f  (T  .)C^^\ 
n-1  n-1  ^  n-1  i_j  p  n-1  n-1 


p=0 

q-1 


(3b) 


C  =  (R  -C  )g(T  )  +y  f  (T  )C^P^ 

n-q+1  n-q  n-q  ®  n-q  p  n-q  n-q 

p=0 


(3q) 


,(j) 


By  substitution,  the  derivatives  at  the  sampling  instants,  can  be  eliminated, 
and  an  expression  of  the  following  form  is  obtained; 

C  .  =  a  (t)R  +a  ,(t)R  ,+.  ..+a  (t)R  +b  (t)C  +.  .  . +b  C  (4a) 

n+1  n  n  n-1  n-1  n-q  n-q  n  n  n-q  n-q 
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where  the  and  are  functions  of  the  past  sampling  interval  durations  T.. 
Example:  Let  G(s)  in  Figure  1  be 


G(s)  = 


K 


s(s+l) 

The  output  during  the  n^^ 


sampling  interval  then  satisfies  the  equation 


c(t)  =  KE 


n 


(t-t  )  - 1  +  e 
n 


+  c  +C 
J  n  n  L 


(t  <t<t  , ) 
n - n+1 

(4b) 


which  corresponds  to  Equation  (la).  Differentiating  once  we  obtain 


c(t)  =  KE 


n 


l_e 

n 


(4c) 


Setting  t  =  in  Equations  (8)  and  (9)  one  obtains 


=  KE 

r  T  -  1  +  e'^^n 

+  C 

+  C 

l-e"'^n 

n 

L  n  J 

n 

n 

- 

C  ^  =  KE  , 
n+1  nL 


-T 

1-e  n 


•  -T 
+  C  e  n 
n 


(4d) 

(4e) 


From  (4d)  and  (4e)  one  can  also  write,  correspondingly, 

n-lL  n-1 


C  =  KE 
n 


l+e''^«-l 

+  c  , 

+  C  , 

l-e'^^n-l 

. 

n- 1 

n-  1 

C  =  KE 
n 


e'Tn-l 

J  n-1 

(4f) 

(4g) 


Terms  containing  and  can  be  eliminated  using  Equations  (4d),  (4f)  and 

(4g).  Then  substituting  E  =  R  -  C  ,  we  obtain: 

^  n  n  n 


-T  1  _e"^'^n+T(i-i) 
K(T  -1+e  «)- 

n 


+e 


-\in‘^^n-l/“I  r,  rp  rr,  . 


(  n-1 


(5) 
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This  is  the  desired  equation,  a  difference  equation  with  time  varying  coeffi¬ 
cients,  which  yields  the  output  at  any  sampling  instant  as  a  function  of  the 
values  of  the  input  and  output  at  the  two  previous  sampling  instants  and  of  the 
lengths  of  the  two  past  sampling  intervals.  Two  past  samples  are  required 
since  this  is  a  second  order  system. 

It  should  be  noted  that  no  restriction  has  been  made,  in  the  development 
of  this  equation,  on  the  manner  of  obtaining  the  time  intervals  T  .  Jury  and 

g 

Mullin  treat  the  case  where  the  sampling  pattern  is  periodic,  repeating  it¬ 
self  every  k  sample,  i.  e.  ,  T^  =  In  tins  case,  one  obtains  difference 

equations  with  periodic  coefficients  which  can  be  solved  by  means  of  z-trans- 
forms. 

g 

Hufnagel  considers  the  use  of  difference  equations  in  the  form  here  de¬ 
veloped  to  perform  a  sample-by-sample  transition  between  two  regions  where 
the  sampling  is  periodic  but  at  different  frequencies.  Standard  z-transform 
techniques  are  used  in  the  first  region  and  the  variable  coefficient  difference 
equations  are  used  to  obtain  required  initial  conditions  to  make  possible  the 
use  of  z-transforms  in  the  second  region. 

If  the  sampling  period  is  signal-dependent,  say  on  the  magnitude  of  the 
error  at  the  preceding  sampling  instant,  the  equations  can  still  be  used  to  ob¬ 
tain  aui  exact  transient  solution  to  this  nonlinear  problem,  proceeding  sample- 
by-sample  in  the  time  domain.  While  this  procedure  is  clearly  very  laborious, 
it  can  be  easily  mechanized  for  solution  on  a  digital  computer. 

3.  ANALYSIS  USING  MODIFIED  z-TRANSFORMS 

5  19 

Modified  z-transform  techniques  have  been  applied  by  several  writers  ' 
to  the  analysis  of  sampled -data  systems  with  a  cyclic  variation  in  sampling 
rate.  It  should  be  noted  that  in  such  systems  the  corresponding  samples  in  each 
cycle  are  separated  by  the  same  time  interval,  as  illustrated  in  Figure  2. 

Since  a  fundamental  periodicity  exists  in  cyclic -rate  systems,  it  is  not  sur¬ 
prising  that  z-transform  techniques  are  applicable.  Two  approaches  to  the 
problem  will  be  illustrated  below. 
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SAMPLING  PATTERN  IN  A  CYCLIC -RATE  SYSTEM 

FIGURE  2 

Consider  first  the  system  of  Figure  3  which  does  not  include  a  hold  cir¬ 
cuit.  Let  the  sampling  times  be  t ^  tg. . .  t^  and  the  cycle  time  T  be  as  shown 
in  Figure  2,  where  there  are  n  samples  in  each  cycle  pattern.  Then,  for  a 
system  with  a  basic  sampling  frequency  ,  the  output  can  be  obtained  by  re¬ 
calling  that  the  convolution  sununation  for  a  linear  sampled-data  system  can 
be  expressed  in  the  form 


00 


ERROR-SAMPLED  SYSTEM  WITHOUT  A  HOLD  CIRCUIT 

FIGURE  3 

Each  of  the  samples  in  a  cycle  is  delayed  by  L.  Therefore,  its  response  at 
time  (rT)  can  be  obtained  by  appropriately  advancing  the  input  function  and  de¬ 
laying  the  weighting  function  g(t),  so  that 

oo 

c.(rT)  =  )  g(kT-t.)e(rT-kT+t.)  (7) 

1  1  X 

k=0 
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E^(z)  =  zE(z,  m)  if  =  mT  (12) 

Similarly,  the  quantities  G..(z)  are  related  to  the  modified  z-transforms  by  the 

relation 

G..(z)  =  zG(z,  m),  ^13) 

J  ^ 

In  order  to  express  the  system  output  in  terms  of  the  input  r(t)  one  can  make 
use  of  the  expression 

Ej^(z)  =  Rj^(z)  -  Cj^(z)  (14) 

Substituting  (14)  into  (10)  one  obtains: 
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(15) 


n 


n 


1  ( 


where  6.,  is  the  Kronecker  delta  function,  i.  e. ,  6..  =  1  for  j  =  k  and  6.,  =  0 

Jk  JK  JK 

for  j  5^  k. 


Hufnagel  suggests  that  Equation  (15)  can  be  solved  conveniently  by  con¬ 
sidering  the  G.,  (z)  to  be  components  of  a  matrix  and  R,  (z)  and  G,  (z)  to  be  com- 
jk  k  K 

ponents  of  vectors,  so  that  Equation  (15)  can  be  written 

[l+G(z)]  {C(z)}  =  [g(z)]  {r(z)}  (16) 


where  ^  is  the  unit  matrix.  The  output  can  be  obtained  by  premultiplying  by 


{C(z)} 


I  +  G(z) 


-  -1  - 


G(z) 


{R(z)} 


(17) 


The  outputs  at  the  sampling  instants,  c(t),  are  obtained  by  evaluation  of  the  in¬ 
verse  z-transform  of  the  Cj^(z)  by  conventional  methods. 


It  should  be  noted  that  when  hold  circuits  are  present,  the  matrix  ele¬ 
ments  G.,  (z)  must  be  modified  since  the  hold  periods  are  of  inequal  length. 
J* 

3.  1  Modified  z-Transforms  from  Switch  Decomposition 


The  result  obtained  above  for  cyclic  rate  sampled-data  systems  can  also 

be  obtained  in  an  easily  visualized  manner  by  following  a  procedure  suggested 
19 

by  Tou.  The  cyclic-rate  sampler  of  Figure  3  can  be  redrawn  as  shown  in 
Figure  4  with  n  constant-rate  samplers,  for  a  system  with  n  pulses  per  cycle. 


EQUIVALENT  DIAGRAM  FOR  CYCLIC -RATE  SAMPLER 

FIGURE  4 
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Each  of  the  equivalent  samplers  in  Figure  4  is  represented  as  preceded  by  a' 
fictitious  adveince  element  and  followed  by  a  fictitious  delay  element,  the  ad¬ 
vances  and  delays  being  equal  to  the  respective  time  interval  between  samples. 
If  the  elements  preceding  the  sampler  are  denoted  by  Dj^(s), 

ti-S 

Dj^(s)  =  e  ^ 

and  the  elements  following  the  sampler  by 
Gj^(s)  =  e  ^^G(s) 

then  the  output  can  be  obtained  by  summing  the  individual  responses, 

k=n 

C(s)  =  ^  Gj^(z)EDj^(z)  (18) 

k=0 

Taking  the  z-transform  of  Equation  (18)  we  obtain 

n-1 

C(z)  =  ^  Gj^(z)EDj^(z)  (19) 

k=0 

The  terms  EDj^(z)  are  obtained  from  the  relation 

E(s)  =  R(s)  -  C(s)  (20) 


by  substituting  C(s)  from  Equation  (18),  multiplying  by  D(s)  and  taking  z-trans- 
forms.  As  a  result  one  obtains  a  matrix  equation 

I^GD(z)  {eD(z)}  =  {rd(z)}  (21) 

where  the  elements  of  the  column  vectors  [ed(z)}  and  {rd(z)}  are 
EDj^(z)  and  RDj^(z),  respectively  and  the  elements  of  the  matrix 

GD(z)  are  of  the  form; 


6..  +  G.D.  (z) 
Jk  J  k 


where 


for  j  =  k  and  6 


jk 


=  0  for  j  5^  k. 


The  solution  follows  lines  analogous  to  those  presented  in  the  previous 
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section  and  produces  the  modified  z-transform  of  the  output,  C(z,  m). 

It  should  be  noted  that  the  application  of  z-transforms  is  limited  to  cy¬ 
clic  variations  in  the  sampling  pattern,  i.  e. ,  there  must  be  a  fundamental 

periodicity  in  the  sampling.  The  technique  has  been  extended  to  the  multiple 

21  9 

sampler  case  and  to  the  systems  with  finite  pulse  width.  If  a  multiple 

20 

sampler  system  is  subjected  to  a  random  input,  it  can  be  shown  that 
"Tuning"  of  one  of  the  samplers  to  operate  at  an  appropriate  cyclic  variable 
rate  can  reduce  the  mean-square  error  below  that  obtainable  with  constant 
rate  sampling. 

4.  SKIPPED-SAMPLE  PROBLEMS:  TRANSFER  MATRIX  FORMULATION 

In  some  systems  samples  may  be  missed  occasionally,  in  an  otherwise 
periodic  constant  sampling  rate  system.  If  the  samples  are  missed  at  random, 
the  problem  may  be  treated  by  considering  the  sampling  frequency  itself  to 
an  appropriate  random  variable.  ^  If  the  samples  are  missed  at  regular  inter¬ 
vals,  the  problem  reduces  to  that  of  a  cyclic -rate  system,  which  can  be  handled 
by  the  methods  introduced  in  previous  sections.  In  this  section  we  shall  con¬ 
sider  an  elegant  solution  to  this  problem,  using  a  matrix  approach  to  find  the 
"skip-sampled  components"  of  a  signal,  as  presented  by  Friedland.  ^ 

Consider  the  sampling  pattern  shown  in  Figure  5,  which  also  shows  one 

skip- sampled  component  of  the  basic  sampled  signal.  The  signal  is  defined 

only  at  the  sampling  instants  nT.  The  component  x.(nT)  of  x(t)  is  given  by  taking 

3 

every  Nth  sample  beginning  with  the  jth  (as  shown  in  Figure  5  for  N  =  3,  j  =  1). 


X(nT) 

0  12  3  4 

Xo(nT) 

0  3  6  9  12  -*-nT 

A  SAMPLED  SIGNAL  AND  THE  ZERO-th  SKIP-SAMPLED 
COMPONENT  IF  N  =  3 

FIGURE  5 
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Thus,. the  N  skip  Sampled  components  are  defined  as 

f  x(nT)  ,  n  =  (j+kN)  (j  =  0, 1, . . .  N-1) 


x.(nT)  = 


n  (j  +kN) 


k  =  0,  1, . . . 


Then  the  signal  x(nT)  is  represented  as 


A 

x(ur)  =  ^  Xj(nT) 
j=0 

The  z -transforms  of  these  components  are  given  by 

00  00 

X^(z)  =  ^  x^(nT)z‘"  =  ^x{(j+kN)T}x'^'^^^^ 


A 


Equation  (25)  may  be  considered  as  a  column  matrix  representing  every  signal 
x(t)  in  the  system  in  terms  of  its  skip- sampled  components,  i.  e. , 

■Xo(*) 

X,(z) 

{X(z)}  =  : 


<26) 

The  input-output  relation  for  an  arbitrary  linear  time-varying  element 
in  the  time  domain  can  be  expressed  as  follows  (at  the  sampling  instants  only): 


■  y(0)  ■ 
y(T) 
y(2T) 


h(0,  0) 
h(T,  0) 
h(2T,  0) 


h(T,  T)  0 

h(2T,  T)  h(2T.  2T) 


x(2T) 


where  h(t,  t)  is  the  weighting  function  for  the  system,  x(t)  the  input  and  y(t) 
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the  output  vectors  respectively,  and  H  =  ^h(nT,  kT)j  is  the  transmission 
matrix  of  the  system  in  the  time  domain.  Now,  since  the  skip- sampling  is 
periodic,  occurring  every  N  samples,  we  have 

h(nT+N,  kT+N)  =  h(nT,  kT) 


The  periodicity  makes  it  possible  to  partition  the  transmission  matrix  into 
(NxN)i  submatrices  and  the  column  vectors  into  (Nx  1)  submatrices.  Then, 
Equation  (27)  is  written  as  follows: 


y(0) 

y(NT) 

y(2NT) 


h(NT) 

h(2NT) 


h(NT) 


x(NT) 

x(2NT) 


In  terms  of  these  submatrices,  the  input-output  relation  can  also  be  written  as 


y(nNT)  =  ^  h(nNT  -  kNT)  x  (kNT) 


which  is  the  matrix  equivalent  of  the  ordinary  convolution  summation  in  the  sta¬ 
tionary  case.  By  defining  the  z-transforms  of  the  matrices  of  Equation  (29), 
and  modifying  them  with  the  appropriate  powers  of  z,  one  obtains  the  matrix 
representation 

{Y(z)}  =  [h(z)]  {x(z)}  (31) 

where  H(z)  is  the  transmission  matrix  in  the  frequency  domain,  or  the  "Trans¬ 
fer  matrix".  For  a  non-time-varying  system,  since  h(nT,  kT)  =  h(nT-kT),  the 


transfer  matrix  becomes: 


Ho(z) 


Hj(z) 


H„(z) 


Hj(z) 


H(z)  = 


«0<^> 


where  Hj(z)  is  the  pulse  transfer  function  corresponding  to  the  appropriate 
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skip-sampled  transmission  matrix.  If  the  component  is  time-varying  but 
has  "zero-memory",  then  H(z)  reduces  to  a  diagonal  matrix,  A  =  diag|^a(0), 
a(T),  a(2T), ...  a  [(N-1)T]  j  .  The  construction  of  the  matrices  is  best  illus- 

4 

trated  by  an  example  from  Friedland’s  paper. 

Consider  the  system  of  Figure  6,  where  the  sampler  closes  for  two 
sampling  instants  and  skips  one  sampling  instant.  Thus,  N  =  3.  The  skip- 


H  6 


SKIP-SAMPLED  SYSTEM  AND  SAMPLING  PATTERN 

FIGURE  6 

sampling  and  hold  can  be  represented  by  a  separate  matrix  which  takes  into 
account  the  fact  that,  in  general,  the  length  of  the  hold  period  will  depend  on 
the  sampling  pattern.  In  this  simple  case,  the  hold  is  described  by 


for  sampling  instants  1,  4,  7, . . .  If  the  integration  is  combined  with  the  plant 

transfer  functiop  — ,  the  hold  can  be  represented  by  the  matrix 
s 


The  plant,  including  the  integration  from  the  hold  circuit,  is 


G^s) 
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so  that 


G(z)  = 


-1 

z 

■ 

G2(z) 

Gj(z) 

= 

Gj(z) 

Go(z) 

G2(z) 

-  G2(-) 

Gj(z) 

Go(z) 

where  the  (z)  are  obtained  by  summing  the  respective  skip-sampled  compo¬ 
nents  0,  1  and  2  of  G(z)  and  are  given  by 


z  ^  -h  2z  ^ 
(l-z"3)^ 


GjCz)  =  ^2li£i 


The  output  of  the  system  can  be  obtained  by  writing  the  overall  transfer  matrix 
for  the  system. 

K(z)  =  I  -  (l  -t-GH(z))  (38) 

where  I  is  the  unit  matrix.  The  response  is  then  evaluated  from 

C(z)  =  K(z)R(z)  (39) 

0 

4.  1  Skipped-Sample  Problems:  Inversion  by  Long  Division 

A  less  elegant  but  conceptually  very  simple  method  of  obtaining  the  output 
of  a  skipped -sample  system  is  based  on  the  well  known  inversion  of  the  output 
C(z)  by  long  division  to  obtain  a  series  in  z  The  coefficients  of  the  series 
are  the  values  of  the  output  time  function  c(nT)  at  the  sampling  instants. 

Consider  the  representation  of  Figure  7,  where  the  skip-sampler  is  re¬ 
presented  by  a  "negative  sample  generator",  which  cancels  pulses  at  the 
"missed"  sampling  instants.  Then,  the  transform  of  the  system  error  can  be 
written  as 


R(z)-HV(z) 
- i  +G(z) 


'  '  1  +  G(z)  (40) 

where  V(z)  is  the  transform  of  the  "cancelling  sampler"  output,  i.  e. , 
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(41) 


V(z)  =  ^  e(nT)z 
n=k 

where  k  represents  the  values  of  the  sampling  instant  n  at  which  samples  are 
omitted.  If  samples  are  omitted  periodically,  then  k  =  Nn,  where  N  repre¬ 
sents  the  particular  skipped  sample;  for  example,  N  =  3  if  every  third  sample 
is  skipped. 


REPRESENTATION  OF  SKIPPED -SAMPLE  PROBLEM 

FIGURE  7 


Now  Equation  (40)  could  be  used  to  obtain  e(nT)  by  long  division,  if  V(z) 
were  known  analytically.  If  it  is  not  known,  the  process  can  still  be  carried 
out  sample  by  sample.  At  each  sampling  instant  where  an  error  sample  is 
omitted,  a  term  corresponding  to  V(z)  is  added  to  the  numerator  (or,  more 
conveniently,  to  the  remainder);  at  all  other  sampling  instants  V(z)  is  zero 
and  therefore  can  be  omitted  entirely. 

If  a  zero-order  hold  circuit  is  present,  the  procedure  must  be  modified 
since  the  output  of  the  hold  circuit  is  not  zero  following  a  "missing"  sample, 
but  rather  will  be  the  same  as  that  due  to  the  last  existing  sample  pulse. 

This  process  can  be  represented  by  effectively  sampling  the  last  existing  value 
again,  i.  e. ,  if  the  n-th  sample  is  omitted,  we  let 

v(nT)  =  -{e[(n-l)T]  -e(nT)}  (42) 

Then 

e(nT)  =  e(n-l)T 


and  the  sampler  output  is  equal  to  the  pulse  applied  at  (n-l)T  once  again.  The 
hold  output  will  then  remain  constant  through  the  "missed  pulse"  period. 
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Therefore,  the  function  v*(t)  can  be  represented  as 


v*(t)  =  ^{a(nT)-e[(n-k^)T]}  6(t-nT) 


(43) 


n 


where  n  represents  the  number  of  the  missed  samples  and  the  number  of 

intervals  since  the  last  existing  sample.  The  application  of  this  method  will 

0 

be  illustrated  with  a  simple  example. 

Let  the  "plant"  G(s)  be  of  second  order, 

G(s)  =  ^ 


s(s+l) 

with  no  hold  circuit  used.  Then 


G(z)  = 


1.  264  z 


-1 


(l-z~l)(l-0.368z“^) 


(44) 


(45) 


Let  the  input  be  a  unit  step  applied  at  t  =  0.  Then 

1 


R(z)  = 


1-z 


-1 


and  the  error  transform  is 

R(z)-V(z) 


E(z)  = 


(rk)- 


V(z) 


1-K3(z) 


1  + 


1.264z’l 


(l-z“^)  (1-0.  3682  ^) 


Equation  (47)  reduces  to 


^  (z  -  0.  368  z)  -  Viz)  (z-1)  (z  -  0.  368) 
(z^  -  0.  104  z  +  0.  368) 


(46) 


(47) 


(48) 


Now,  let  the  sampling  pattern  be  such  that  every  third  sample  is  missed,  as 

_2 

shown  in  Figure  5.  Then  the  first  term  in  V(z)  is  e(2T)z  and  V(z)  can  be 
omitted  for  the  computation  of  the  first  two  samples,  as  shown  below 

1  -0.  264z"^ 

2 


z'^-O.  104z  +  0.  368  'z2-  0.  368  z  [+  f(z)V(z)l 

*■ -  J 


0.  104z  +0.368 


-  0.  264z  -0.368 
-0.264Z  +0.027-0.097  z"^ 


-  0.395+0,  097  z 


-1 
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Now,  since  there  is  no  sample  at  the  instant  (2T).  we  add  v(2T)  to  the  remain¬ 
der  and  continue.  Starting  from  the  last  remainder,  we  add  -0.  395  z^(z-l) 
(z-0.  368)]  : 

0  z~^  -  0.  443  z~^ 

z^-0.  104  z  +  0.  368  I  T  0.  395  +  0.  097  z“  J 

-  0.  395  +  0.  540  z~^  -  0,  146  z~^ 

-  0.  443  z'l  +0.  146  z"2 

-  0.  443  z~^  +0.  046  z~^  -  0.  163  z~^ 

+  0.  100  z”^  +0.  163  z”^ 

Thus  the  first  four  terms  of  e(nT)  are  then: 

e(0)  =  1.0,  e(T)  =  -0.264.  e(2T)  =  -0.  395,  e(3T)  =  -0.  443 

(49) 

The  output  is  then  computed  from  the  relation 

c(nT)  =  r(nT)  -  e(nT)  (50) 

It  should  be  noted  that  while  this  method  does  not  give  the  output  in 
closed  form  directly,  it  is  not  limited  to  periodically  skipped  samples.  The 
missed  samples  can  be  designated  by  an  arbitrary  function  of  (nT)  in  advance, 
or  they  can  be  computed  as  a  function  of  system  performance  at  the  preceding 
intervals.  Thus,  for  example,  the  sampler  may  be  designed  in  such  a  manner 
that  it  samples  at  the  instant  (n+l)T  only  if  e(nT)  >  €,  i.  e. ,  if  the  error  ex¬ 
ceeds  certain  limits.  Similarly,  the  sampler  operation  may  be  governed  by  a 
performance  criterion  such  as  the  sum  of  a  number  of  terms  of  the  mean- 
square  error  sequence.  However,  since  the  solution  proceeds  sample-by¬ 
sample  from  t  =  0,  it  is  probably  limited  to  the  investigation  of  the  first  few 
sampling  intervals  if  it  is  not  to  become  excessively  time  consuming. 

5.  THE  STATE  VARIABLE  FORMULATION  OF  SAMPLING  PROBLEMS 

A  formulation  of  sampled  data  problems  using  the  concepts  of  state  vari¬ 
ables  and  state  transition  matrices  has  been  presented  by  Kalman  and  Bertram, 
The  method  is  primarily  applicable  to  systems  with  periodic  sampling  (either 
constant- rate,  cyclically  varying  rate,  or  multi-rate  systems),  and  has  been 
applied  to  the  analysis  of  systems  with  random  sampling.  ^  ^  As  with  other 


10 
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methods,  the  generalization  to  arbitrary  sampling  periods  makes  the  method 
extremely  laborious.  In  fact,  Kalman  and  Bertram  emphasize  that  their  ap¬ 
proach  is  designed  for  use  with  digital  computers. 

The  formulation  of  the  method  begins  by  dividing  the  components  of 
linear  sampled  data  systems  into  three  groups:  (a)  continuous  dynamic  ele¬ 
ments  (characterized  by  linear  differential  equations);  (b)  discrete  dynamic 
elements  (characterized  by  linear  difference  equations);  and  (c)  sample-and- 

r 

hold  elements  which  represent  the  transitioii  from  discrete  to  continuous 
information. 


The  state  of  the  continuous  elements  is  represented  by  the  (yx  1)  matrix 
[x^(t)]  where  7  is  the  sum  of  the  orders  of  the  differential  equations  governing 
these  elements;  the  elements  of  [x^(i)]  the  values  of  the  continuous  quan¬ 
tities  x^  and  their  derivatives  at  time  t.  The  state  of  the  discrete  dynamic 
elements  at  time  t  is  the  set  of  numbers  required  at  the  start  of  a  cycle  of 
computation  to  solve  the  particular  set  of  difference  equations.  They  are  de¬ 
noted  by  a  6 -vector  x^(t)]  there  are  5  such  numbers  required. 

The  state  of  the  sample  and  hold  devices  is  denoted  by  a  column  matrix 
whose  elements  are  the  values  of  the  outputs  of  these  devices  at  time  t.  This 
is  denoted  by  [  x®(t)l  and  is  assumed  to  be  of  rank  (ctx  1). 


The  state  of  the  entire  system  is  then  obtained  by  combining  the  three 
state  vectors  into  a  single  (nx  1)  matrix  which  the  authors  call  [  x  (t)]  .  i.  e. , 


{x  (t)}  = 


x''(t) 

x'^(t) 

x®{t) 


(51) 


The  transitions  from  samplirig  instant  to  sampling  instant  are  then  writ¬ 
ten  by  considering  the  state  transitions  of  each  type  of  element  to  be  a  matrix 
acting  on  the  state  variables. 


The  derivation  of  the  transition  matrices  is  beyond  the  scope  of  this  re¬ 
port,  but  the  method  will  be  illustrated  with  one  of  Kalman  and  Bertram's 
examples.  Consider  the  system  of  Figure  8,  which  includes  discrete 
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compensation  as  well  as  a  continuous  "plant". 


ILLUSTRATIVE  PROBLEM  FROM  KALMAN  AND  BERTRAM 

FIGURE  8 


The  state  variable  vector  for  the  continuous  element  is  given  by 


(52) 


If  the  discrete  compensator  is  described  by  the  difference  equation 


(53) 


then  a  new  variable  x  0(^1,)  can  be  defined  by 

u  iC 


(54) 


Thus,  the  state  of  the  element  at  time  tj^,  given  the  input  e(ij^)  is  clearly 
described  by 


_X°(tk)]  =[x3\.l)_ 


(55) 


The  state  of  the  sample -and -hold  element  is  simply  its  output,  which,  at  the 
sampling  instants,  equals  its  input: 
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[xV]  '[*4<v]  -[«v]  (56) 

In  the  above  equations,  t^^  is  the  kth  sampling  instant  after  t  =  0,  and 

"  \+l  “  V* 

Now  it  can  be  shown  that  the  state  of  the  complete  system  at  time  t. 

K  •  X 

is  described  by  the  equation 


{x  (tk+l)}  =  {x  r(tj^ 


where 


{x  (ti  4 )}  is  the  state  vector  at  time  t,  , 

'•  k+1  ^  k+1 

is  the  system  transition  equation  from  instant  tj^  to  instant 

^k+1 

r(tj^)  is  the  input  at  time  tj^  (a  scalar) 

y(Tk)  J  is  a  vector  including  the  effect  of  the  input  on  the  state 

transition. 

For  the  particular  system  of  this  example,  the  transition  matrix  is  given 


)i- 


:Tkx,_-Tk 


.  ■'(I'k)  ]  ■ 


.o(l-e-Tk) 
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The  top  two  rows  of  the  matrices  in  Equations  (58)  and  (59)  govern  the  state- 
transitions  of  the  continuous  elements;  the  third  row  concerns  the  discrete 
compensator  and  the  fourth  row  the  zero-order  hold. 

While  the  formulation  of  the  equations  in  matrix  form  is  very  elegant, 
it  should  be  noted  that  Equation  (57)  is  merely  a  concise  statement  of  the 
difference  equation  method  outlined  in  a  previous  section  of  this  report.  The 
matrix  formulation  has  particular  appeal  in  the  stationary  case.  Even  in  the 
periodically  varying  sampling  rate  case,  when  the  problem  is  stated  in  terms 
of  difference  equations  with  periodic  coefficients,  the  "transition  matrix"  be¬ 
comes  stationary  if  attention  is  focused  on  the  transition  corresponding  to  a 
complete  period  of  sampling  operations. 

6.  OTHER  METHODS 

6.  1  Continuous  Approximations 


If  the  sampling  frequency  u  is  sufficiently  high  compaured  to  the  band- 

s 

pass  of  the  plant  G(ju),  then  the  sample-and-hold  operation  can  be  approxi¬ 
mated  with  one  or  two  terms  of  its  frequency  domain  representation.  A  zero- 
order  hold  is  described  by 


T 

-ju)  Y 

H(ju)  =  Te 


sinu  T/2 
tjT/2 


(60) 


Therefore,  the  pulse  transfer  function  of  a  sampler  followed  by  a  zero-order- 
hold  and  plant  G(ju)  is  given  by 


.  -jT(u+nUg)/2  sin  (u+nWg)T/2 

HG*(ju)  =  -  )  Te  - ^ ^  G(ju  +j 

T  ^  (w+nu„)T/2 

o 


n=-oo 


nu)  ) 
s 


Since  for  low  frequencies  ^  higher  harmonics  are  assumed 

negligible,  this  transfer  function  reduces  to 

HG*(ju)  =  e'^^'^^^Gtiu)  (62) 

This  approximation  is  the  foundation  of  the  most  common  approximation  meth- 
2  12  13 

ods.  Linvill  '  has  suggested  the  extension  of  the  approximation  by  using 


ASD  TDR  62-35 


23 


one  or  two  of  the  additional  sidebands  due  to  sampling.  Clearly,  the  diffi¬ 
culty  with  applying  this  type  of  approximation  to  the  variable  sampling  rate 
problem  is  that  one  obtains,  at  the  very  least,  continuous  control  systems 

with  variable  time  delays.  The  solution  of  problems  of  this  type  is  also  at 

T  16 

best  approximate.  Other  continuous  design  methods  '  suffer  from  similar 

17 

limitations.  Rubtsov  has  presented  a  careful  analysis  of  the  validity  of 
some  of  the  approximation  methods. 

6.  2  Slowly-Varying  Sampling  Rate  Systems 

If  the  sampling  rate  varies  slowly  in  comparison  to  the  dominant  sys¬ 
tem  time  constants,  other  approaches  are  possible.  (The  concept  of  "slow 
variation"  can  be  related  to  the  change  in  the  sampling  rate  during  the  time 

required  for  the  impulse  response  of  the  system  to  decay  below  a  specified 

18 

value. )  Tartakovskii  has  applied  the  methods  used  by  Zadeh  for  the  study 
of  stability  of  time  varying  linear  equations  to  the  stability  of  difference  equa¬ 
tions  with  varying  coefficients.  This  method  makes  it  possible  to  apply  z-trans- 
form  techniques  to  obtain  approximate  stability  conditions  for  slowly  varying 
systems. 

If  the  sampling  rate  varies  slowly  it  may  also  be  considered  piecewise 
constant.  Ordinary  z-transforms  can  be  used  for  the  solution  during  a  con¬ 
stant-rate  period,  and  the  difference  equations  of  the  system  solved  step-by- 

0 

Step  at  transitions  of  sampling  frequency  to  establish  initial  conditions. 

7.  SUMMARY  AND  CONCLUSION 

This  report  has  presented  a  survey  of  methods  presently  available  for 
the  analysis  of  sampled  data  systems  with  a  variable  sampling  rate.  A  num¬ 
ber  of  exact  methods  can  be  used  for  systems  where  the  sampling  rate  varies 
periodically,  since  the  inherent  periodicity  makes  the  use  of  z-transforms 
possible. 

The  general  time-dependent  sampling  interval  problem  can  be  solved 
by  writing  the  state  transition  equations  of  the  system  from  interval  to  interval. 
These  equations  can  be  expressed  in  concise  form  using  matrices,  but 
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except  in  extremely  simple  cases,  they  are  too  complex  and  time  consuming 
for  solution  except  with  a  digital  computer.  A  number  of  the  methods  pre¬ 
sented  are  practical  for  evaluation  of  system  response  during  the  first  few 
sampling  instants,  even  if  the  sampling  intervals  are  given  by  arbitrary  func¬ 
tions  of  time. 

\Vhen  the  sampling  rate  is  dependent  on  system  performance  the  prob¬ 
lem  becomes  nonlinear.  Recursion  equations  for  solution  of  the  nonlinear 
difference  equations  can  still  be  written  and  the  above  remarks  apply.  Most 
of  the  other  methods  surveyed  in  this  report  cannot  be  extended  directly  to  , 
the  nonlinear  case. 

While  stability  evaluation  was  not  discussed  in  the  report,  Kalman  and 
Bertram^ ^  indicate  the  extension  of  stability  criteria  from  constant-rate  to 
periodically  varying  sampling  rate  systems.  No  general  results  appear  to 
be  available  for  the  arbitrary  variable-rate  sampled-data  system. 
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Unclassified  Report 

This  report  presents  several  techniques 
which  can  be  used  for  the  analysis  of 
sampled-data  systems  with  a  non-constant 
sampling  period.  It  is  shown  that  the  appli¬ 
cation  of  z-transform  techniques  is  limited 
to  cases  where  the  sampling  pattern  repeats 
periodically.  Several  special  case^,^ 
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including  the  cyclically-varying  sampling 
period  and  the  "Skipped  sample"  problem 
are  outlined.  The  most  general  methods 
available  are  based  on  a  direct  solution  of 
the  system  difference  equations.  These  equa¬ 
tions,  while  time  consuming  to  solve,  do 
make  possible  the  evaluation  of  transient 
response  sample-by-sample.  The  advantages, 
limitations  and  possible  extensions  of  the 
various  methods  are  outlined.  The  report 
includes  a  number  of  simple  examples  and 
an  extensive  bibliography. 
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